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Abstract 

A Monte Carlo generator has been constructed to simulate the reaction e^e~ 
7 + 47r, where the photon is assumed to be observed in the detector. Only initial 
state radiation is considered. Additional collinear photon radiation has been in- 
corporated with the technique of structure functions. Predictions are presented for 
cms energies of IGeV, 3GeV and lOGeV, corresponding to the energies of DA$NE, 
BEBC and of i?-meson factories. The event rates are sufficiently high to allow for a 
precise measurement of R{Q^) in the region of Q between approximately IGeV and 
2.5GeV. For the construction of the program we employ isospin relations between 
the amplitudes governing r decays into four pions and electron positron annihilation 
into four pions. Estimates of the kinematic breaking of these isospin relations as a 
consequence of the vr" ~ vr'' mass difference are given. 



1 Introduction. 



The precise determination of the cross section for electron positron annihilation into 
hadrons over a large energy range is one of the important tasks of current particle physics. 
The results are relevant for the analysis of electroweak precision measurements which are 
affected by the running of the electromagnetic coupling from the Thompson limit up to 
Mz- Also the interpretation of the anomalous magnetic moment of the muon depends 
critically on these data. Last not least the measurement of the energy dependence of R(s) 
is one of the gold plated tests of QCD and allows for a precise determination of the strong 
coupling constant. 

Depending on the energy region different techniques for the measurement of R{s) 
have been applied up to date. At low energies, say from the two pion threshold up to 
roughly two GeV, exclusive channels are collected separately, for higher energies inclusive 
measurements start to become dominant. For energies below m,- isospin invariance and 
CVC have traditionally been used to predict r decays from electron positron annihilation 
[0, H, 1^, 1^ . Clearly this strategy can be inverted pending irreducible uncertainties from 
isospin violation and radiative corrections . At high energies a multitude of final states 
is present and only inclusive measurements have been performed. 

To cover a large range of energies, results from many different experiments and colliders 
have to be combined, and energy scans have to be performed to obtain the full energy 
dependence. An attractive alternative is provided by the upcoming $— and B— meson 
factories which operate at large luminosities, albeit at fixed energies. Events with radiated 
tagged photons give access to a measurement of R over the full range of energies, from 
threshold up to the CMS energy of the collider. For events with tagged photons the 
invariant mass of the recoiling hadronic system is fixed by the photon energy which 
provides an important kinematic constraint. The usage of photons observed at extremely 
small angle with respect to the beam has been investigated in |0, . In this case final state 
radiation as background is practically irrelevant. However, in practice photon detectors do 
not cover this region. At large angles a careful study of initial versus final state radiation 
has to be performed. 

To arrive at reliable predictions including angular and energy cutoffs as employed by 
realistic experiments, a Monte Carlo generator is indispensable. For hadronic states with 
invariant masses below two or even three GeV it is desirable to simulate the individual 
exclusive channels with two, three up to six mesons i.e. pions, kaons, etas etc. which 
requires a fairly detailed parametrization of the various form factors. 

In principle initial and final state radiation would be required for the complete sim- 
ulation. Such a program has been constructed for the two pion case 0. There it is 
demonstrated that suitably chosen configurations, namely those with hard photons at 
small angles relative to the beam and well separated from the pions, are dominated by 
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initial state radiation. In fact, this separation is possible |]T0[ even when operating the 
(f) factory DA$NE on top of the resonance where direct radiative decays cannot be 
ignored. 

In the present paper we continue this project with the construction of a generator for 
the radiative production of the four pion final state, including the uo{-^ 37r)7r channel. 
This mode contributes a large fraction of the rate with invariant masses between one and 
two GeV. The energy region between 1.5 GeV and 2.5 GeV is difficult to access directly 
with current electron positron colliders. At the same time the experimental uncertainties 
are relatively large. This motivates the special effort devoted to this range. 

The Monte Carlo program discussed in the present paper is specifically devoted to 
cms energies up to approximately 10 GeV. It is constructed in a modular form such that 
the parametrization of the hadronic matrix element can easily be replaced by a more 
elaborate version. Different final states with three, four or five pions or kaons can be 
included. The present parametrization of the hadronic matrix element follows closely the 
form suggested in correcting only some minor deficiencies. The four pion amplitude 
is assumed to be dominated by p/ — > vrai plus a direct coupling pi —> pun and exhibits 
the proper behavior in the chiral limit. 

The plan of this paper is as follows: In the next section the formalism for the decom- 
position of the differential cross section into a leptonic and a hadronic tensor is presented. 
Results for partially integrated distributions are recalled which can easily be used to ar- 
rive at simple estimates for the rates. In section 3 isospin relations are derived between 
the amplitudes for four pion production from a virtual photon and those accessible in r 
decay. The relations between these four matrix elements contain the well known iden- 
tities between the corresponding rates and provide in addition important constraints on 
the differential distributions. In section 4 the influence of the 7r° — vr^ mass difference 
on the relations obtained in section 3 is discussed. The ingredients of the ansatz for the 
hadronic amplitude together with a comparison between the model prediction and the 
data for a variety of distributions are presented in section 5, the complete description of 
the hadronic amplitude with all the model parameters is collected in the appendix. A 
description of the Monte Carlo generator is given in section 6, together with a few char- 
acteristic distributions. In particular we investigate the influence of angular cuts and the 
addition of collinear radiation. Section 7 contains a brief summary and the conclusions. 



2 The radiative return. 

Hard photons observed at small angles relative to the electron or positron beam and at 
the same time well separated from charged particles in the flnal state can be used to 
reduce the effective center of mass energy at electron positron colliders. Performing a 
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detailed analysis of the angular and energy distributions for the 77r"'"7r~ final state it has 
been shown that initial and final state radiation can be reasonably well separated . 
For the four pion case we therefore restrict the discussion to initial state radiation only. 
The matrix element for the production of an arbitrary hadronic final state corresponding 
to the diagrams Fig.|l| are given by 




Figure 1: Diagrams contributing to the process e^e ^ 7 + hadrons (only initial state 
radiation included). 

The matrix element of the hadronic current 

jr = jr(gi,-,gn) =<%i),...,Mgn)iJr(o)io> (2) 

has to be parameterized by form factors to be discussed below. For the two pion case the 
current 

JT''^ = - 9;) F,^ (g2) (3) 

is determined by only one function, the pion form factor F2Tr. 
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For the three pion case the matrix element of the electromagnetic current is restricted 
by current conservation and negative parity to the form 



and is dominated by the cu resonance. The matrix element for the four pion case will 
be discussed in sections 3 and 4. 

The differential rate can be cast into a product of a leptonic and a hadronic tensor 
and a corresponding separation of the phase space 



da^- L^^H''- d^2{p+ +P-; Q, k) d^n{Q; qi,...,qn)^ , (5) 

OS ZTT 

where d^n{Q', Qi, ■ ■ ■ ,Qn) denotes the n body phase space with all statistical factors in- 
cluded. 

The leptonic tensor L^j, is process independent, the modeling of hadronic physics 
enters the tensor Hjj^i, = J^"^{J^"^)* only. The leptonic tensor is symmetric, and hence it 
is only the real symmetric part of H^i, which enters. It has the following form 
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{kp_){kp+) 



g'"' [p-p+ + 



g5 1 + P+P+ + P-P- 



1 



(6) 



The coUinear and soft photon singularities proportional to 1/(1 ± cos 6*^) and 1/ are 
evident from these expressions, as well as the l/Q^ enhancement for small Q^. Integrating 
the hadronic tensor if^;^ over the hadronic phase space one gets 



j J^rUrr d^r.{Q-.qi,...,qn) = -^(Q,Q.-9,.Q') R{Q') (7) 

where R{Q^) is (7(e+e'~ — > hadrons) / apoint- 

The additional integration over the photon angles (the azimuthal angle is integrated 
over the full range and the polar angle within 9min < < T^ — dmin ) leads to the differential 
distribution 



2 da /(s2 + g4) 1 + cos^^in {s - Q^) \ 

which can be used to calculate the event rate observed for realistic photon energy and 
angular cuts (see Tab.l). 
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Event rates 


Collider 




Integrated luminosity, fb ^ 


■ — ^° 


a -TO 

"min ' 


^min 10 


DA$NE 


1.02 


1 


13 • 10^ 


12 ■ 10^ 


10 ■ 10^ 


B — factory 


10.6 


100 


4- 10^ 


3.5 ■ 10^ 


3 ■ 10^ 


B — factory 


10.6 


100 


2.7- 10^ 


2.3 ■ 10^ 


2.0 ■ 10^ 



Table 1: Estimated number of radiative events e~^e~ — >■ hadrons + 7 for different center 
of mass energies from Eq.(H). In the first two rows by hadrons we mean vr+vr" and the 
minimal photon energy is 0.1 GeV. The third row is obtained for a continuum contribution 
in the region 2 GeV < y/Q^ < 3.7 GeV assuming a constant R = 2.4. 



3 Isospin relations. 

The emphasis of this paper is towards hadronic final states consisting of four pions and a 
photon. Before entering a discussion of a model dependent parametrization of the form 
factors (see section 5) we recall the constraints from isospin invariance. They relate the 
amplitudes of the e^e^ 27r"'"27r~ and e"'"e~ — >■ 7r"'"7r~27r'^ processes and those for r 
decays into 7r^37r'' and 7i^27i^7i^. The amplitude of the r decay into an arbitrary number 
of hadrons plus a neutrino is given by 



Mr 



Gp 
V2 



cos 9c V (p^) 7" (1 - 75) M (pr) 



(9) 



with 



4 = Ja i(lu-,1n) = <h{qi),...,h{q„)\J^{0)\0> 
and (0) = d "iaU ■ (10) 

We use the same letter J for the operator and its matrix element and restrict our consid- 
erations to the Cabbibo allowed vector part of the hadronic current. 

It leads to the differential distribution 



dV 

dQ^ 



;ii) 



with 



I J^J-* d^„iQ;q,,...,q^) = — {Q,Q,-g,,Q') W{Q^) . (12) 

Note the relative factor of 2 between the definitions in Eq.(|^) and Eq.(0). 

Ignoring the issues of isospin breaking and radiative corrections, the electromagnetic 
current can be decomposed into an isospin singlet piece and a part transforming like the 
third component of an isospin triplet: 



r'^ = ^ + ^ (13) 

72 3^2 



whereas the charged current generating r decays is given by 



J- = i= (Ji-zj2) . (14) 

Final states with an even number of pions are produced through the isospin one part 
only, whence 



V2 J7(27r) = J^-(27r) (15) 

and R{Q'^) = R^{Q^) for two pion final states. 

A similar relation for the four pion final state is easily obtained as follows: the tran- 
sition from the vacuum to four pions is mediated through a current (0) of the form 



J^i (0) = ^ J d'^qid^q2d'^q3d'^q4 Jf, (tti ■ vr2) (tts x n^) , (16) 

where vTj = n{qi) denotes the pion field in the Cartesian basis. The letter J is again 
used both for the operator and the function in the integrand which corresponds essen- 
tially to the transition amplitude. The function = J^(gi, 52, Q's, ^'4) is symmetric (anti- 
symmetric) with respect to the interchange of qi and q2 {qs and q^). The combination of 
pion fields relevant for the transition to four pions with total charge -1 and respectively 
is given by 

1 r 

J'f, (0) = 4 y d^qid*q2d^q3d^q4 ■ (nln^ + 7rf vr^ + tt^tt^) (^VTg tt^ - n^n^) 
4(0) = \ Jd%d%d%d% J^-{nin, +7,^71^ + 7,y,){n^n,- .(17) 
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Taking the matrix element of these operators between vacuum and the states 
(7r^(]9~'')7r~(p~)7r'^(pi)7r°(p2)| etc. one immediately arrives at 

(TT+TT-TrOTrOl^JIO) = MPuP2,P^,P~) 
(7ri+7r^7rf7r^|jJ|0) = J^^{pt ,P2 ,Pi ,Pi) + JM ,P2 ^Pt ^Pi) 

+ Jli{pt^Pl^Pl^P2) + Jt^{Pl,Pl,P2^P2) 

(7r"7r$*7r^7r^|J^ |0) = J^(p2,P3,P",Pi) + ^/.(pi, Ps, P", P2) + Jf^{Pi,P2,P~ ,P3) 
(7rf7r^7r+7rVu |0) = J^{p^,P2,PuP°) + Jf^ip^ ,PuP2,P^) , 



which connects r decay and electron positron annihilation. It is clear from Eq.(|T8|) that 
only one matrix element, namely the one for (+ — 0), needs to be programmed and 
the remaining ones can be obtained by relabeling arguments. Interference terms between 
the two partitions |jT2|, |13|, ^ (3,1) and (2,1,1) are present in the differential distributions. 



For the integrated rates induced by the currents one obtains 



R{+-0 0) = ; +0)=A + ^B 

R{+ + - ~) = A + B ■ i?^ (- 0) = i (A + 5) , (19) 

with 



2-71 f 

^ = ~'q2 J ^^(?i'92,g3,g4) (gi,g2,g3,g4) c?$4(Q;gi,...,g4) 

^ = ^^('^''*^^i'^2,g3,g4)^;i(gi,g3,g2,g4)) rfl>4(Q;gi,...,g4) , (20) 

consistent with the familiar relations between r decays and e+e^ annihilation into four 
pions: 



0) = ^ i?(+ + - -) 
+0) = i i? (+ + - -) + i? (+ - 0) . (21) 
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4 The TT^-TT^ mass difference and the isospin relations. 

All the relations obtained in the previous section are strictly applicable only in case all 
pions in the final states have the same mass, which is obviously not true. The relatively 
large (about 3.6%) vr^-vr^ mass difference will affect the -R(Q^) ^ K^i^Q"^) relation even if 
the relations Eq. (|15|) and Eq. (|18|) still hold. The CVC hypothesis and the assumption that 
transitions to an even number of pions in the final state are described by the iso-vector 
current are well established experimentally |T^. It is thus natural to assume that the 



relations between currents hold and the question is up to what precision we can ignore 
the 7r^-7r° mass difference. This should give at least an indication of the size of these 
"kinematic" isospin violations. We will address this issue here first for two pion states in 
some detail and subsequently indicate the size of possible effects for the four pion case. 

In order to estimate the size of kinematical isospin violations we first start from the 
assumption that the cross section cr(e"'"e~ —* vt+tt^) is well measured and the squared 
form factor is extracted through 



\P(Q")\" = - vr+7r-)(g^) , (22) 

= (1 — 4m^_/(5^)^''^) and used to predict the r~ — >• Ut-tt'tt^ decay rate (We ignore the 
electroweak correction factor S^w — 1.019.). The size of the corrections depends critically 
on the details of the assumptions. 

If the form factor and the form of the current 

J7(27r) = V2 j;(27r) = (g^,, - g^,,) (23) 
remain unchanged, the integral for the r rate is given by 



r(r Ut-ti 71^^ 

COS^ Or 



8ml 



dQ' f 1 - ^ 
V mi 



/3!o f 1 + + 3 (3.0 



m. 



(m^ — mp)^ 



\Fm\,m 



with 



(m_+mo)^\ / {m_ — moY 
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1/2 



(25) 



The second term is due to an S-wave contribution and can be eliminated by replacing 
Eq.(^ by 



V2 J^-(27r) = (^gi,, - g^,, - ^'q^ ^'(^i - ^^)) ^(^') • (26) 

Numerically the contribution to the integral of this latter term is tiny - nevertheless we 
shall adopt the second choice. 

This purely kinematic modification^ raises the prediction for the r decay rate by 0.86%. 
It seems, however, plausible, that the energy dependent width of th p- meson, which is 
present in the form factor, has to be modified accordingly]^ leading to an effective increase 
of T p by 0.74%. The two effects nearly compensate in the integral. Hence the relation 
between r~ z/^7r~7r° partial decay width and the e^e~ —>■ 7r"'"7r~ cross section based on 



T(t- 



3 cos^ 9^ 




(Tee 



TT^TT ) 



(27) 



would only be corrected by 0.06%. 

However, as shown in Fig.^, a sizable dependence of the ratio of the two spectral 
functions (V^Q'^) ~ ) is expected, with a reduction approximately 0.74 % 

close to the peak of the p resonance and enhancements at the tails. 

At present, however, r data provide an important input for the prediction of the QED 
coupling at the scale of Mz and the hadronic contribution to g-2 |TB 



Let us, for 

the moment, assume that the aforementioned kinematic effects are indeed present. If the 
e"'"e~ cross section is deduced from r data through Eq.(pS|) 



1 dr{T Ut-TT tt 

3 cos^ e, 



Q' 1-^ U + cr(e+e- TT+TT-) (28) 




-'^This is the strategy adopted in the generator TAUOLA jlj which was written when data were not 
precise enough to require the incorporation of mass corrections in the form factor. 
^J.K. thanks M. Davier for drawing his attention to this point. 
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1.05 




Vs [GeV] 

Figure 2: The ratio of the two spectral functions (see text for explanation). 





only (3^ modified 


and Tp modified 


0.28 < E [GeV] < 0.81 








1.0163 


1.0088 




1.0116 


1.0027 


0.32 < E [GeV] < 1.777 






Aa{Mz) 


1.0130 


1.0058 


1.0096 


1.0016 



Table 2: Kinematic correction factors for the predictions of and Aa from r data. 
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and the contributions to g-2 and a{Mz) are evaluated without kinematic corrections of 
phase space and form factor the former are overestimated by 0.58% and 0.16% respectively 
(Tab.2). 

The situation is more complicated for the four pion case. All Air modes have different 
numbers of 7r°, whence the phase space is different for each of the mode. Moreover, for 
a comparison between R{Q^) and R^iQ^) one has to integrate over four particle phase 
space and it is impossible to obtain a simple analytical result. To estimate a size of the 

effect we integrated the quantities R{ hO) etc. according to Eq.(|ll]) using the current 

discussed in the next section once assuming that all masses are equal to m_ and once 
taking the real masses. For the mass corrections of the integrals we find 

(-000) : 5.0% ; (--+0) : 2.4% ; (+ + ) : 0% ; (+-00) : 4.6% (29) 

Mass effects alone thus modify the integrated version of the Eq.(|ll|) to 



) 

where 

r(-OOO) = T{t- z/^7r"37r°) ; r(- - +0) = r(r" z/^27r"7r+7r°) 

r(+ + -) = 2 / (i - «;) + 2 «;) + -, 

r(+ - 00) = 2 r,^^ / f 1 - ^ V f 1 + Ri+ - oo) dQ' . (3i) 

In principle this correction depends on the form of the current which will be specified 
in the next section. However, since ai and u dominance and the qualitative behavior of 
the matrix element is well established this result will be applicable also to other forms 
of the current. Additional effects could arise from modifications of the widths in Breit- 
Wigner functions. The corrections are thus just indicative and detailed studies would be 
required which are beyond the scope of this work. 

It is also instructive to consider the mass effects on the differential rate for the 
e+e~ — > 27r°7r+7r~ process. In Fig.| we plot the ratio TZd = ""^"""^ , where cxi ((J2) is 



1 

1.050 



1.024 



-000) 

-+0) 



2r(+ + 
^r(+ + 
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0.12 — 



O.l - 



0.08 



0.06 



0.04 



0.02 




Figure 3: The ratio 71^ 



0-2 



where cxi ((T2) is e+e 27r°7r 7r+ cross section 



calculated for true pion masses (with all masses equal to m_). 



the e^e" 27r'^7r^7r+ cross section calculated for true pion masses (with all masses equal 
to m_). The difference amounts again up to a few percent. The relations between the 
differential r decay rates and the e~^e~ cross sections obtained in the previous section 
will be violated at that level of accuracy even if Eq.(0) holds. However, to test these 
predictions experimentally more precise measurements of the e~^e~ — An cross section are 
required, where now a typical systematic error amounts to roughly 15 %. 



5 The hadronic current. 

As one can see from section 3 it is enough to construct only the hadronic current for 
the (H — 00) mode, while the other ones can be obtained using the relations Eq.( 



Its construction was based on |]TT| with some small changes allowing for preserving the 
relations Eq.(|l^). The basic building block of the current contains a part built on the 
assumption of Oi vector dominance plus an u exchange contribution. However only by 
adding an /o contribution one can recover the proper chiral limit |l^. The complete 
current can be written as a sum of these three contributions 



^pO-.27TO-K + 7T- - Tfli + ^fo + ■ (32) 

They are depicted schematically in Fig.^ and described in detail in the Appendix. Here 
we present some numerical tests of the current and a comparison between results obtained 
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p .p 




P (0 / 7U 
->■ ^ — > f — >- 



7U 




Figure 4: Diagrams contributing to the hadronic current. 



using the current Eq.(^) and experimental data. One should add that the parameters of 
the model are kept as in [O] even if in principle they should be re-fitted as the current 



is a bit different from the original one of |]TT| and new and improved data have became 
available. This, however, is beyond the scope of this paper. 

Starting from tests of the code of the current, first one can check if the Monte Carlo 
program reproduces the known ||T^ analytical results of the partial r decay widths in the 
chiral limit 



71 



r (r ^ z/^27r-7r+7r°) cos^ 



rrir 



1 /1009 



15 \2TTf^l 128 V 96 



TT 



(33) 



73 



r (r ^ UrTT Stt^ 



cos^^c ( Y 1 f 2 _ 437\ 
15 Utt/^J 256 V "Is"/ 



(34) 



Eq.(P5D differs from Eq.(38) of |T^. It was not discovered there that the analytical 
result is wrong as the tests were done at 0.2 % (one sigma) precision level and the difference 
amounts to 0.34 % . Tau decay rates can be obtained in two different ways. One way 
is just their direct calculation. The second one is by using the known relations to the 
e^e~ — ^ 27r+27r^ and e~^e~ 7r+7r^27r° cross sections [§, |lll. The results of the numerical 
tests are summarized in Tab. 3, where CVC means the decay width was obtained through 
its relation to the e'^e~ cross sections. As one can see the results of the test performed 
at 0.02 % accuracy level are quite satisfactory. This agreement gives confidence in the 
numerical stability of our program. 
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Decay mode 


MC result (direct) 


MC result (CVC) 


analytical result 


71 


0.026786(4) 


0.026786(4) 


0.026788 


73 


0.015998(3) 


0.015996(2) 


0.015999 



Table 3: Comparison between analytical and Monte Carlo results in the chiral limit. 



Mode 




present model 


experiment 


Br {r- z/^27r-7r+7r0) 


3.11% 


4.33% 


4.20(8)% 


Br (r^ UrTT'uj (vr^vr+Tr^)) 


1.20% 


1.48% 


1.73(11)% 


Br (r- Ur7i-37i^) 


0.98% 


1.14% 


1.08(10)% 



Table 4: Branching ratios of r decay modes. Results of |Tl[ and the present current are 
compared to experimental data. 



Now we can test the physical predictions of the current. Let us start with r decay 
branching ratios. The results are summarized in Tab. 4, where we put for completeness 
also the results presented in [|TT|] and the experimental results []TB[. The agreement of the 
predictions of the current Eq. (|32D with the experimental data is satisfactory for r~ 
i/t-t^'Stt^ decay mode. Comparing however the results for the r~ — >■ i't-2tt^7t^tt^ and 
UriT^uj (tt^tt^tt'^) modes it seems that the u part of the current does not represent 
the data well, even if the total branching ratio for the v^2'k^'k^'k^ decay mode 

agrees with the data. Again the results of the partial decay widths were obtained as in 
the chiral limit by direct calculation and checked by relating them to the simulated e^e~ 
cross sections. Agreement was found within statistical errors proving that the code of the 
current fulfills the CVC relations Eq.(^) (if integrals are performed with m_ = mo). 

In the next step the predictions for the e"'"e~ 27r+27r^ (Fig.|^) and e+e" 7r+7r^27r'' 
(Fig.^ cross sections are compared with data. The plots taken from |19| contain data sets 
from OLYA |2g, ND ^ CMD [2^, SND |2|], CMD2 [|g plus results of the Orsay 

The data have a typical systematic error of 



25, 26, and Frascati 



groups. 

about 15% (shown in plots only for some of the data sets) and we have thus decided not 
to refit the parameters entering the current as the agreement between the Monte Carlo 



and the data is acceptable. The modification of the current of |11] we performed leads 
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Figure 5: Comparison of data (left figure, see text) and predictions (right figure) obtained 
using the current Eq.(p2D (filled squares) and those of |Tl| (empty squares) respectively 
for cT(e"'"e~ 27r+27r^). 




Figure 6: Comparison of data (left figure, see text) and predictions (right figure) obtained 
using the current Eq.(P^ (filled squares) and those of |TT| (empty squares) respectively 
for aie'^e" 7r+7r^27r^). 
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to significantly better agreement between the theoretical predictions and the data in the 
mode e^e'^ 7r+7r~27r°. Considering the agreement between the r decay rates and the 
prediction we conclude that the model reproduces the data well even if the description of 
the uj part of the current is not completely satisfactory. 



6 The Monte Carlo program.^ 

The idea behind the structure of the Monte Carlo program is to allow for a simple addition 
of new final state modes into the program and for a simple replacement of the current (s) 
of the existing modes. The program thus exhibits a modular structure. For the generation 
of the four momenta of the mesons no sophisticated method of a variance reduction was 
applied. This slows down the generation, but could be accounted for if a faster Monte 
Carlo generator would be required. It has, however, the advantage of being universal 
and no change of the variance reduction method is required with each modification of 
the hadronic current. The process to be simulated by the program in its final stage is 
e"'"e~ — i> 7 + hadrons with an exclusive description of final states, even if till now only 
TT+TT", 27r'''7r+7r~ and 27r"'"27r~ hadronic final states are implemented. The LL radiative 
QED corrections were taken into account using structure function method as developed 
in ||3l| and limited to the initial emission only. In fact the program can run in one of 
two modes (chosen by a user) one with coUinear radiation and one without it. Hard large 
angle photon emission is limited to initial state radiation, which is justified by where it 
was demonstrated for the vr+vr" hadronic state that the contribution from the final state 
emission as well as the initial-final state interference can be reduced to a negligible level 
by applying suitable cuts. 

The generation of the multi particle phase space is based on the following representa- 
tion of a Lorentz invariant phase space 



dLips„+i (P, /c, gi, g^) = ,^ dQ^ ... dQ.^_i dLipsa (P, Qi, k) dLipsg (Qi, Q2, Qi) 

(27r) 

... dLipS2 (Q„-2,Qn-l,gn-2) dLipSa ((5n-l,gn-l,gn) , (35) 

where P is a total four momentum of the photon + hadrons state (it is not equal to 
the sum of the initial e^e~ four momenta as we allow for an additional initial coUinear 
emission), k is the photon four momentum, qi ... g„ are the four momenta of the hadrons. 



Qi = qi+ ... + Qn and Qi = Qi^i - qi-i for i = 2, ... n-l , (36) 



^The program is available upon request from the authors. 
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and dLips2 (^1,^2, ^3) is a two body phase space 



dLips2 {ki, k2, ks 



1 xHkikiki) 



(37) 



327r2 A;? 



with A (a, b, c) = + 6^ + — 2ab — 2ac — 26c and dQ^ the ka sohd angle. 

The generation flow is the following: first the collinear radiation is generated and a 
four momentum P of the visible final state is calculated and then boosted to its rest 
frame (RF) (a CM of the visible final state). This part is omitted in the mode running 
without collinear emission and then P is a sum of the electron and positron four momenta. 
In the second step the visible hard photon four momentum is generated. Its energy is 
generated flat even if the energy distribution is govern by a competition between soft ~ 
and a complicated resonant spectrum which depends on details of the hadronic current. 
Of course this way the generator is not extremely efficient, but it is more universal and 
will work equally well with all reasonable modifications of the hadronic current. The 
photon polar angle is generated in P RF with the distribution accounting for collinear 
emission peaks, which are the same for all hadronic modes as the initial state is ever the 
same. Its azimuthal angle again is generated with a fiat distribution in P RF. Then a 
chain of a generations of , i = 2, ...,n — 1 follows (Qi is fixed when we generate the 
photon four momentum). They are generated fiat within their allowed limits 



where Q\ is a zeroth component of the Qi four momentum. At the end a generation 
of the solid angles fli, i = 1, ...,n — 1 follows. They are generated fiat ( fli in Qi RF). 
All generated four momenta are then transformed with the use of the proper boosts and 
rotations into the CM system of the initial e~^e~ particles. The distribution of the events 
~ |A^P, where is a matrix element of a given process is obtained by means of the 
rejection method. The cross section is calculated using Eq.(|]) both for weighted and 
unweighted event samples. The unweighted events are stored in a file when requested. 

Let us now discuss the cuts, which reduce the contribution of the final state emission 
to the cross section to a negligible level allowing thus for extraction of R{Q^) from the 
measurement of the e~^e~ hadrons + 7 cross section. We recall that in it was 
shown that the following set of angular cuts 




(38) 



cutsl : 



(7° <9^< 20' 



1° 



or 160° <9^< 173°) and 30° < 9^ < 150° , (39) 
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( 9^ {6t^ ) is the photon (pion) polar angle) fulfils this requirement for the e~^e~ tt~^7[~'j 
cross section. It reduces, however, the observed cross section significantly. This starts to 
become dramatic, when one runs at energies well above 1 GeV. The following set of cuts 



cuts2 : (7° <9^< 20° and 30° < 9^ < 173" 

or (160° <9^< 173° and 7° < 9^ < 150°; 



(40) 



also reduces the contribution from final state radiation to a negligible level due to the fact 
that the pions and photon are well separated as in the previous case. 
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Figure 7: The differential e'^e~ Att^ cross sections at beam energy 1.5 GeV with 
minimal photon energy equal 0.1 GeV with no cuts on pions angles and 7° < 9^ < 173° 
( no cuts) and two sets of angular cuts Eq.(p9D (with cutsl) and Eq. (|40|) (with cuts2), 
where q = \fQ^ is an invariant mass of the 47r system. 



At the same time the cross section reduction is much smaller, especially for higher 
beam energies and higher energies of the observed photons. The effect of the two sets 
of cuts on the cross sections with four pions in the final state is presented in Fig.|^ for 
2Ebeam = 3 GcV. For higher beam energies the effect of the cross section reduction is 
much higher and at 2 Ebeam = 10 GeV the cross sections for cuts specified in Eq. ([39|) is 
reduced almost to zero. For the cuts specified in Eq.(|y) the results are presented in Fig.|| 
and the reduction remains tolerable. 

From Fig.^one concludes that one can measure R{Q^) at a B factory in the interesting 
region of \/Q^ between 1 GeV and 2.5 GeV through measuring the e^e~ — >■ Airj cross 
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Figure 8: The differential e^e~ — > 4717 cross sections at beam energy 5 GeV with minimal 
photon energy equal 0.2 GeV with no cuts on pions angles and 7° < 9^ < 173° ( no cuts) 
and angular cuts of Eq.(^OD (with cuts2), where q = \fQ^ is an invariant mass of the 47r 
system. 







Event rates 




Integrated luminosity, fb ^ 


27r+27r-7 


iT^^-K^-n 7 


1 GeV 


1 


1.04- 10^ 


1.13 ■ 10^ 


3 GeV 


1 


4.66 • 10^ 


5.72 ■ 10^ 


10 GeV 


100 


1.86 • 10^ 


2.33 ■ 10^ 



Table 5: Estimated number of radiative events e+e^ ^ 47r + 7 for different center of 
mass energies. The minimal photon energy is: 0.05 GeV (first row), 0.1 GeV (second 
row), 0.2 GeV (third row). The angular cuts of Eq. (|40D were applied. 
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section using Eq.(^. This measurement should have an an accuracy much better then 15 
% , which is now the typical experimental error in that region, allowing for a reduction 
of the error in the calculation of the photon vacuum polarization. From Table 5 it is clear 
that the error would be dominated by systematics and not by statistics. 

One may even restrict photon and pions detection angles to the central region, e.g. 
25° < < 155° and 30° < Oj^ < 150° respectively if a minimal angle of 20° between 
photon and charged and neutral pions is required in order to suppress final state radiation 
and to clearly separate neutral pions and the photon. With this cut one obtains {y/s = 
10 GeV and C = 100 /fe"^ ) a rate of 1.17 ■ 10^ events with 27r+27r~7 and 1.53 ■ 10^ events 
with 27r°7r"'"7r~7. 




Figure 9: The effect of the collinear radiation on the differential e~^e~ An^y cross sections 
at beam energy 0.5 GeV (left picture) and 5 GeV (right picture) with minimal photon 
energy equal to 0.05 and 0.2 GeV correspondingly. In the mode with collinear radiation 
minimal invariant masses of the 47r + 7 systems of 0.95 GeV and 9.5 GeV were required. 
q = \fQ^ denotes the invariant mass of the 47r system. 

Additional collinear emission always present in the real experiment reduces slightly the 
cross sections as shown in Figs.^ for two different modes and two different beam energies. 
Its actual size depends on the cuts on the invariant mass of the 47r + 7 system. The 
effect is similar for different energies and for both charge modes. 
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7 Summary. 



A precise value of the cross section for hadron production in electron positron annihilation 
at low energies is one of the important ingredients for a reliable prediction of the anomalous 
magnetic of the muon and the electromagnetic coupling at high energies. As an alternative 
to a direct measurement at the relevant energy one may use initial state radiation to 
reduce the effective energy of electron positron colliders, exploiting the large luminosity 
of "factories" and accessing thus a continuum of hadronic final states. 

With this motivation a Monte Carlo generator has been constructed to simulate the 
reaction e^e~ — 7 + 47r, where the photon is assumed to be observed in the detector. 
The hadronic matrix element has been taken from [|ll|]. Isospin relations between the 
amplitudes governing r decays into four pions and electron positron annihilation into four 
pions have been found which allow to determine all four modes after the amplitude for 
the 7r+7r~27r° channel has been fixed. The kinematic breaking of these isospin relations 
as a consequence of the vr" - 7r° mass difference has also been investigated. 

The program is constructed in analogy to the one @ simulating e"^e^ ^ 7 + 27r. 
However, it does not include final state radiation from the charged pions. Additional 
coUinear photon radiation has been incorporated with the technique of structure functions. 
Predictions are presented for cms energies of IGeV, 3GeV and lOGeV, corresponding to 
the energies of DA<I>NE, BEBC and of i?-meson factories. 

Even after applying realistic cuts the event rates are sufficiently high to allow for 
a precise measurement of RiQ"^) in the region of Q between approximately IGeV and 
2.5GeV. 

The model predictions are compared to recent data from electron positron colliders. 
Once more accurate data become available, the modular structure of the program will 
allow for modification of replacement of the hadronic current in a simple way. 

Acknowledgments. The work by H.C. was supported by BMBF-POL-239-96, the 
work by J.K. by BMBF under grant BMBF-057KA92P. H.C. is grateful for the support 
and the kind hospitality of the Institut fiir Theoretische Teilchenphysik of the Karlsruhe 
University, where this work was carried out. 

Appendix: The description of the hadronic current. 

In this appendix we give a complete definition of the hadronic current used in this 
paper. Recalling Eq. (|3^) 

rpO-+27r"7r+7r- = ^ai + ^'fo + , (41) 
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we will define here its ingredients. 

Denoting the pions four momenta as follows: q'i(vr°), 52(71"°), ^3(7?") and (3'4(vr"'") one 
gets the following contribution from the part containing an ai exchange 



r^i (gi,g2,g3,g4) = 

(93,52, 91,54) +f;^^ (93,91,52, 54) - f^i(94,92,9l,93) - f^i(94,9l,92,93)- (42) 

The structure of this current follows from the form of the Lagrangian of the pan 
a ■ [n X p) ) and piTTT (~ p- (tti x 712) ) interactions and the function is defined as 

H 

f^:, (91, 92, 93, 94) = Ci T,(g2) r'^-(Q,g-gi) 

BWa,{{Q-qiy) T,x{Q-qi,Q-qi-q2) 

BW,{iq3 + q,f) rt(g3 - 94), (43) 



where we have followed the notation of [|ri[]; Q = 9i + 92 + 93 + 94- 

Assuming that 9i = 9I = 93 = 9I = m"^ i.e. that pions have equal masses one finds 

f^i(9i,92,93,94) = 

Ci fpnTT ((93 + 94)^) faipAQ"^! {Q - 9i)^) faipAiQ - 9l)^ (93 + 94)^) 

T,(g2) BWa,{iQ-qif) BW,{iqs + q,y 



92(93-94) ^^/^(9i + 92)(93 - 94) , (Q9i)(92(93 -94))^ 



(93 - 94)" + 9r 7A V2 - Q 



,(44) 



(Q-9i)^ " V Q\Q-qi? 

where /p^^r and fa^pir describe pTrvr and aipvr vertices appropriately while Tp, BWa^ 



and -BVFp are p, oi and again p propagators. The different choice of the authors of ||TT 
for two p propagators comes from the fact that p couples in a different way in the two 
cases and the 'propagators' absorb in fact some parts of the couplings. The normalization 
constant Ci = was chosen to give the proper chiral limit of the current. The form 
factors were chosen to be constant 

[/p.. ((93 + 94)') /a,p.(g^ (g - 91)') fa,p.m - 91)', (93 + 94)')] = o.ss, 

and the validity of that assumption was proved [|n[] in a limited range of 
(1.1 GeV < < 2.2 GeV ). 
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The contribution to the current from /o exchange reads 



r/o(?i,92,g3,g4) 



C2 fppfo (Q^ (?3 + g4)^ (gi + ^2)^) fpnTT ((ga + ^4)^) f/onn {{qi + q2)' 

Tp{Q') Tp[{q, + q,f) BW {{q, + q,)' 



iqs-q.r-Q^^^^ 



(45) 



where /pp/,, and ffonn describe pp/o and /oTttt vertices appropriately, while BWf^ is an 
/o propagator. 

Again (72 = —^75^ follows from the chiral limit and the form factors are kept constant 

fppfo {Q\ (QS + 14?, (QI + 12?) fpn. ((ga + Qif) ffonn + 52)') = 0.38 . 

The contribution coming from so called anomalous (containing u exchange) part of 
the current reads 



ri;(gi, 92, 93,^4) = ^ 1475.98 GeV"^ 12.924 GeV^^ 0.266 ml 
v2 

[gfFi(gi, g2, gs, ^4) + g2^i(g2, gi, gs, q^) + g3^2(gi, 92, gs, 54) - g4^2(gi, 92, g4, gs)] , (46) 

where 

Fiiqi, q2, qs, q^) = BWpM\ iQ - ^2)') [{qsiQ - ^2)) (g2g4) - {q^{Q - ^2)) (g2g3)] (47) 
and 

^2(gi, g2, g3, 94) = SWp,^(g', (g - 92)') [(g4(g - 92)) (gig2) - (gi(g - 92)) (g2g4)] 

+ BW,M\ - qi?) [(?4(g - qi)) (gig2) - (g2(g - gi)) (gig4)] (48) 



and g^^ = 1.55. It is changed from its original value g^^ = 1.4 in [|Tl| to reproduce the r 
decay rates r~ ^ z/^27r~7r+7r° and r~ — z/^7r~u; (vr" 71+71*^). 

For completeness we list here all propagators lll| required for the current: 



Tp {Q') 
with 



BWsiQ^ ^p, Tp) + PiBWsiQ^ m^,, T^J + /32 W3(Q^ m^,, J 

l + /5l + /52 



(49) 
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and the numerical values set to 



with 



= 0.14 GeV 
rup = 0.773 GeV = 0.145 GeV 
rup^ = 1.35 GeV Tp^ = 0.3 GeV Pi = 0.08 
rup^ = 1.7 GeV Tp^ = 0.235 GeV p2 = -0.0075 ; 



EWp{Q') - BWs{Q\mp,Tp)+pBWs{Q\mp,,tp,) ^ 



rfip^ = 1.37 GeV = 0.510 GeV P = -0.145 ; 



where 



m: 



ai 



9{Q^) ' 



and 



Q 39 fi5 

^(Q^) = 1.623 + 10.38 - ^ + ^ for > (m„, + 

y(g') = 4.1 _ 1 _ 3.3 (g2 _ j ^ ^ g (^q2 _ 9^2 

for g^ < (mai + m^)^ 

= 1.251 GeV r„j = 0.599 GeV ; 
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with 



m/o = 1.3 GeV = 0.6 GeV ; 



(58) 



- inipTp 



+ (T- 



m 



p' 



- inipiTp,^ 



1 



6 ((2.2 GeVf - Q^) + 6 (Q^ - (2.2 GeV)' 



:2.2 GeV) 



2\ 2- 



(59) 



with 



m 



pi 



1.7 GeV r 



pi 



0.26 GeV a 



-0.1 



0.782 GeV T,, = 0.0085 GeV . 



(60) 



An additional suppression compared to |TT| above = (2.2 GeV)^ was introduced. 
This modification prevents the unphysical growth of the cross section for very large Q^, 
which originates from the momentum dependent couplings in Eg. (^6]) -Eg. (^8]) . 
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